A note on the Guerra and Talagrand theorems for Mean Field 
Spin Glasses: the simple case of spherical models 

Silvio Franz and Francesca Tria 1 

1 The Abdus Salam International Center for Theoretical Physics 
Strada Costiera 11, P.O. Box 563, 1-34100 Trieste (Italy) 
E-mail: franz@ictp.trieste.it, ftria@ictp.trieste.it 
(Dated: February 2, 2008) 

Abstract 

The aim of this paper is to discuss the main ideas of the Talagrand proof of the Parisi Ansatz 
for the free-energy of Mean Field Spin Glasses with a physicist's approach. We consider the case 
of the spherical p-spin model, which has the following advantages: 1) the Parisi Ansatz takes the 
simple "one step replica symmetry breaking form", 2) the replica free-energy as a function of the 
order parameters is simple enough to allow for numerical maximization with arbitrary precision. 
We present the essential ideas of the proof, we stress its connections with the theory of effective 
potentials for glassy systems, and we reduce the technically more difficult part of the Talagrand's 
analysis to an explicit evaluation of the solution of a variational problem. 
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I. INTRODUCTION 



The mathematical analysis of the low temperature mean field spin glass phase has seen 
enormous progresses in the last few years, after more than twenty years since the proposal 
of the Parisi ansatz that led to its physical understanding |2[. 

The first important progress was achieved by Guerra and Toninelli ||, who, through 
an interpolation among systems of different size provided a temperature independent proof 
of the existence of the thermodynamic limit for the free-energy and other thermodynamic 
quantities in a wide class of mean field spin glass systems. 

The power of the interpolating method was fully appreciated by Guerra [ffl, who com- 
paring interacting mean field systems with suitable paramagnets, was able to write the 
free-energy as a sum of the Parisi expression plus a remainder term, which is manifestly pos- 
itive in the SK model as well as in p-spin models for even p. The method was subsequently 
generalized to deal with diluted spin glass models H, 0] and more recently to analyze the 
Kac limit of finite dimensional spin glasses ||. 

The last, but definitely not the least, step towards the proof of the Parisi free energy 
solution has been performed by Talagrand who, through a highly non trivial generalization 
of interpolation techniques to coupled replica systems, could finally prove the vanishing of 
the remainder in the Guerra formula |9|, [10], |TTf . 

Unfortunately the elegance of the Talagrand approach is somewhat obscured by the math- 
ematical necessity of dealing with many technical details in the course of the proof. 

In this note we would like to present the Talagrand theorem using a physicist's approach, 
emphasizing the main ideas and connecting it with the theory of the glassy effective poten- 



tials discussed in the physical literature |L2], |13|, [L4|, [L5|, |16[| . 

We discuss the case of the spherical p-spin model, where the saddle point equations take 
a particularly simple form and where it is well known that the best Parisi solution is of 
the simple "one step replica symmetry breaking kind" (1RSB). This allows to reduce the 
technically difficult part of the proof to the explicit solution of variational equations that we 
analyze numerically. 

Our paper comprises six sections. In the second section we define the model and we 
set the basic definitions that we use in the rest of the paper. In the third section we 
illustrate the backbone of the Talagrand theorem, specializing it to our simple case. In the 
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fourth section we derive the replica expressions and evaluate numerically suitable "effective 
potential functions" whose positivity implies the validity of the Parisi Ansatz. We conclude 
the argument in section M, where we show that the replica expressions for the potentials 
provide a lower bound for their exact values. We finally draw our conclusions. 



II. SOME DEFINITIONS AND THE GUERRA FORMULA. 

In the spirit of discussing the results in the simplest non trivial example, we focus on the 
case of the spherical p-spin model. 

The model considers N real spin variables <jj, i — 1, ...,N, subject to a global spherical 
constraint Yli=i of = N, and interacting through a p-body Hamiltonian: 

1,7V 

H[<t\ = - Yl J n-iv a h ■ ' ' a i P ■ W 

U<«2<--*p 

The couplings Ji v ..i are i.i.d. Gaussian random variables with zero mean and variances: 



E, 



J 2 



2NP- 



(2) 



and we are not considering the case of a non zero external magnetic field. 

The best Parisi solution is known to be of 1RSB kind with q = 0. The Guerra interpo- 
lating Hamiltonian J| to the level of 1RSB can be written as: 

TV 

H t [a]=VtH[a}-VT^tJ2(hi + h d , i )(T i , (3) 

i=l 

where the interpolating parameter t runs in the interval [0, 1] and hi and hd,i are Gaussian 
i.i.d. fields with zero mean and variances that can be written in terms of a parameter 
q E [0, 1] as: 

E h [hj\ = \<T X E hd [hi] = |(1 - (T 1 ) • (4) 
One then defines the generalized free-energy per spin: 

ft = -Jjj- Ej m Eh [E hd Z t {K h d ,i)] m , (5) 

where Z t indicates the partition function at inverse temperature (3 computed with the Hamil- 
tonian H t , while m is a number in the interval (0, 1]. Notice that f\ is the free-energy of 
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the original model, while fo is the free-energy of a simple paramagnet which can be read- 
ily computed expressing the spherical constraint through a Lagrange multiplier and using 
standard saddle point method. This gives: 

V° = ~\ 5£ {a - m (a - f (i - (i - m)^)) - (a - ^(i - r». 

(6) 

Formula (|5|) makes a contact between the interpolating method and the cavity method. In 
fact it is well known that (|5]) can also be written in the form []TT , 18|, [l9j (and for a recent 
review see also l2C 



/t = E J,f 7 ln 5^ w i E d z t{h ia , h i)d )) , (7) 

7=1 

where w 7 = „^ are random weights derived from "free-energies" / 7 that realize a 

Poisson point process with exponential density p(f) = /Sme^" 1 ? [p5fl , and the h in are drawn 



independently for each i and 7 from the Gaussian distributions defined by (^). The fields 
/ij jC i, hi receive then naturally the interpretation of cavity random fields, whose variances 
have to be determined self-consistently in order to estimate the free-energy f\. 

By simple integration by part, it is possible to see that the free-energy f t can be written 
as i: 

/* = fo + ^=^(1 " (1 " m)q P ) + £ ds TZ S . (8) 

The first two terms together in (|8]) provide the 1RSB replica free-energy / t rep for the in- 
terpolating model. TZ t is the t derivative of a remainder term which can be written as: 

TZ t = £ (1 - m) (q(a, rf - pq(a, r)<r 1 + (p - 1)?)^ + (q(a, 7-)% , (9) 
where q(a, r) denotes the overlap function: 

N 

X 



1 N 



8=1 



The two averages (-)i t and (-)o,t are defined introducing two kinds of interpolating Hamilto- 



nians for two copies of the system with the same interactions [|10||, of the form respectively: 



A? 



H^[a,r] = VtH[a] + VtH[r] - VT^Yl + + + ' ( n ) 



i=l 
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with fields hi and h r id , r = 1,2, distributed according to (f|) with h\ d and h\ d independent, 
and: 

N 

Hl%,r] = sTtH[o] + V~tH[r] - VT=* £ ({h\ + + (h\ + /»» 4 )r,) , (12) 



i=l 



where now, all the fields ft,* and ft|, and hf d are drawn independently with the distribution 



To enlighten the notations it is useful to define averaged partition functions: 



Zl°\h\h 2 ) = E h>l [Zr(h l ,h z ,h d ,!>, 
zf\h) = E W [z^(h,h d ,h 2 d 

and, for a given function of two spin configurations k(a,r), averaged Boltzmann averages: 

E h ij* /_ + °° Va Vrk(a, r) exp (-/3tf t (1) [a, r] ) 



(13) 



a>x(k(a, r)) 



£j (fc(o-,r)) 



E h>l ^ Prfcfo j exp {-(5Hj 0) [a, r] ) 
Z?\h\h*) 



(14) 



(15) 



A? 



Do- = daS IJ2 °i - N 



i=l 



(16) 



With these definitions, we can finally write: 

E h z[ l \h) 2 u x {h{a,T)) 



(k(a,r)) u = 



(fc(<7,r))„ t EE 



— 

^ (1) (^) 2 

^^^(^^^"^(fc^r)) 



(17) 



(18) 



E hlM zt\h\V) m 

It is apparent from the convexity of the function q p for even p that both averages in 
the remainder (^|) are non-negative 0. Talagrand theorem consists in proving that in the 
thermodynamic limit, for an appropriate temperature-dependent choice of the parameters 
q, m, both terms in @ are indeed equal to zero at all temperatures. For t = 1 the free- 
energy fx is independent of q and m. The minimization of the remainder J dt !Z t is therefore 
equivalent to the maximization of the 1RSB free-energy. We will need to consider the t- 
dependent free-energy, for the values of q and m that achieve this maximization. 



III. TALAGRAND THEOREM. 



The Talagrand's analysis of the remainder can be divided in three logical steps: 

1. From the Taylor formula one simply observes that: 

q(a, rf - pq(a, r)q^ 1 + (p - l)g* < ^-^(g(<r, r) - qf , (19) 
and it is evident for p > 2 that: 

q(a,T)> <q(a,r) 2 ; (20) 
to bound the remainder it is then enough to bound the averages ({q{cr, t) — q) 2 ) 1 1 and 

2. In order to get estimates of the averages ((q(a,r) — q) 2 ) lt and (q(a,r) 2 ) 0t , one con- 
siders the t-dependent probabilities Pj >t (qc) = (8(<l(o r , T ) — Qc))j,t for j = 0,1, that 
the overlap takes a value qc, respectively in the ensembles defined by (-) lt and (-) ot . 
The main part of the theorem consists in proving that, given e > 0, for each t < 1, a 
t-dependent positive constant C{t) exists such that, if {qc — q) 2 > C(t)(f t — fl ep ) + e 
and q 2 c > C(t)(f t — f^ ep ) + e, then it exists a positive constant C'(t,e) such that 
Pi,Mc) < exp(-C'(t,e)N) and P , t (qc) < exp(-C"(t, e)N), where f t ep indicates the 
t-dependent free energy calculated with the Parisi ansatz. 



We will explicitely see in section [IV B| that the constant C'(t, e) turns out to be a 



monotonically decreasing function of t, behaving as e(l — t) 2 for t close to 1. 

3. The estimates in the previous point imply, in the large iV limit and fixing any to < 1, 
the existence of a positive t -dependent constant K(t Q ) for which we can write the 
differential inequality: 

d{ft ~/ r) <K(t )(f t -fr), (21) 

valid for every t < t . Since f t =o = fl=oi © proves f t < f t ev for every t < to, which 
by continuity can be extended to t = 1. 

While the first and third points can be easily understood, the proof of the second point 
requires several additional steps. 



A. Analysis of the probabilities Pjj(qc)- 



To estimate the probabilities Pjj(qc) it is useful to relate them to suitable "potential 
functions" for coupled replicas, similar to the ones used in the physical literature to study 
off-equilibrium configurations of glassy systems J[2], [L3], [14], [15], |l6f . The potential functions 
can be defined as the difference between the generalized free-energies of the interpolating 
system for the two replica with and without an additional constraint in the mutual overlap: 
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ViMc) = 



mN/3 



EjlnE h [Z?\h)uMq(<r,T) -q c ))) ' ~ Ej In E h Z ( t x> (h) 



(22 



EjhxE hlM [Z?\h\h 2 )uJ Q {5{q{<j,T)-q c ))) - EjhxE h , M Z\ u \h\h 2 ) 



mN/3 

Notice, that, despite the dependence does not appear in the notation, the Boltzmann aver- 
ages u(-) depend in (|22"| , |2~3"D on the same fields as the adjacent partition functions Z^ and 
Z^ respectively. 

One can then use the property (0) and write: 

J2y=i w^ZKh^ujjiSiqia, r) - q c )) 



c 



T 

— E Jf In 



c) 



E 



w^ZKhJujiSjqjcT, t) - q c )) 



(23) 



(24) 



Again, in formulae (|2l3| , |2~4|) it should be understood that the fields appearing implicitly in 



the uj(-) are the same of the adjacent partition functions Z{. In the weights u> 7 = |=r 



-Pf-y 



-Bf-y > 



the variables / 7 are chosen with exponential densities respectively px(f) = e~f in the case 
labeled by 1, po(/) — e/3m ^ i n the case 0. 

One can now immediately see that Pj it (qc) and Vj )t (qc) are respectively the expectation 
value and a properly normalized expectation of the logarithm of the same random variable 
Xj t t(qc), where the averages are performed over the random couplings and the free-energies 

u 



v^qc) 

P h Mc) = EjjX j>t (qc 



-^EjjJnX^qc) 



(25) 
(26) 



It is possible to show, using techniques explained for example in pTJ, that 1/N hi(X j tt (qc)) 



is a self-averaging quantity and the following "concentration of measure property" holds: 
Prob [X jtt (q c ) > e - N ^Aia)+Ne 1 ^ + Proh [x jjt (q c ) < e -W.*(ffo)-^i] < 2e~^ (27) 

for all e\ > 0. Noting that Xj t (q) is a probability density, integrating it in q on a small 
interval [g^, qc + Ag] and passing to the limit Aq — * 0, from a direct computation one 
obtains: 

PjMc) < 4e"^ + e- Npv *^-* N + e -^Aio)+e,N (2g) 

We thus see that Pj,t{Qc) is bounded by an exponentially small number in N whenever 
Vj,t(Qc) is strictly positive in the limit N — > oo. 

The remaining part of the theorem is then devoted to showing the positivity of the two 
potentials Vjj{qc) f° r every t < 1. To this end one would like to have a lower bound 
provided by the replica expressions [p6| , that are the quantities that we are able to evaluate. 
Unfortunately, with the interpolation method we cannot directly prove the desired inequality 
VjAlc) > Vj^ilc)) where Vjj P (qc), j = 1,0, indicate the replica expression of the two 
potentials. One can however show that the replica expression gives a lower bound for the 
generalized free energies fjjiqc) — Vj,t(Qc) + 2/t- To infer from this bound the lower bound 
for the potentials Talagrand found an ingenious shortcut: one first prove the existence for all 
t < 1 of a positive constant C(t) such that, for all qc ^ q one has V{J p (qc) > ^ or 1° 

such that (qc — q) 2 > C(t)(f t — f[ ep ) one can then write the following series of inequalities: 

ViMc) > V7{qc) + 2fr(qc) ~ 2f t > Vg*{q c ) ~ > • (29) 

An analogous series of inequalities can be written for Vo t t(qc) for all qc ^ where one finds 
VoTilc) > % for some C(t). 

In the next section we use the standard replica approach to evaluate the two potentials 

v;; p { qc ). 

IV. THE POTENTIAL FUNCTIONS IN THE REPLICA APPROACH. 

In order to get quick heuristic estimates of V\ and Vo we resort to the replica analysis, 
postponing to section [V] the discussion on how the expressions we find provide free-energy 
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lower bounds. The readers not interested to the details of the replica derivation can jump 
directly to section [IV B . 

Before entering the discussion of the potential functions, we find it useful to illustrate 
the analysis of the interpolating model for a single copy of the system through the replica 
method. This can be done starting from expression (^), considering m as an integer and 
substituting the log with an n/m-th power (supposed to be an integer): 

ft = -JJf^ EjhiE h [E hd Z t {hi,h d ,dr - -j^^Ej[E h [E hd Z t (h,h d )) m }™ . (30) 

Expanding the powers one sees that we are in presence of n replicas a a of the original 
systems divided in n/m groups, each with m replicas. Together with the spin configurations 
the fields h are also replicated. Replicas a and b in the same group have identical fields 
h" = h\, while replicas a, h in different groups have statistically independent fields hi and 
hfi. The Gaussian distribution of the fields is summarized in the covariance matrix: 



p/2 a=b 
S p/2 q p ~ l a,b in the same group of m replicas (31) 
a,b in different groups 

At this point the calculation follows the usual rails of replica analysis of mean-field spin 
glasses 0). Introducing the overlap matrix Q ab one finds: 



fr = - iim Up. { if E + i V M E H *Q* + ^ ln det Q } 



(32) 



where S.P. denotes saddle point evaluation over the parameters Q ab . The saddle point 
equations for Q read: 

^fQlb 1 + ^yVa* + ^(Q-'U = 0. (33) 



The values of H ab that make the saddle point equations ( p3|) independent of t, are such that: 

H ab = \Q p - h x ; (34) 
P%Q P ^ = -(Q- 1 )* ■ (35) 

This is the choice needed to minimize the remainder (ft). 



The further analysis of expression (^) is standard |22| and will not be reproduced here. 
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We pass now to the slightly more involved expressions for the potentials Vj tt (qc)- We give 
full details for the potential Vx t t(qc)- An analogous procedure can be followed for Vo t t(qc)- 
The expression for the "constrained free energy" fxj(qc) is given from (|22| ) by: 

2 



c) 



mN(3~ 



Ej\nE h {Z^(h)coMq^r) - q c )) 



1,2 a 



,n (36) 



nN/3 



\nE 



J,h ra ,h r d a 



+oo 



n Va * ex p 



We now consider two copies, indexed by r = 1,2, replicated n times. The fields h r,a and h s,b 
are equal if a and b belong to the same group of ~ replica, while they are independent in 
the other case for all choices of r and s; the fields h r d a are independent for each replica ra. 
In this case we define a field covariance matrix TC ra ,sb = E[h™h!l ] + E Wi "d^t dl wri ich can be 
conveniently recast in the form H a ^ = Ti ra ,rb (t = 1,2), A Qj b = Ti Ta , s b (r ^ s — 1,2). The 
form of these matrices is given by: 

p/2 a=b 
H a b = { p/2 q p ~ l a,b in the same group of m/2 replica (37) 
a,b in different groups 

p/2 qP- 1 a=b 
A a b = \ p/2 q v ~ x a,b in the same group of m/2 replica (38) 
a,b in different groups 

As for the single replica system, the best choice for the values of the parameters q and m 
that appear in the fields' variances is the one that maximizes the free energy of the "free" two 
replica system, that is, the interpolating system without the additional constraint q±2 = qc- 



As explained in []12[ the replica computation of the constrained free-energy is analogous 
to the unconstrained case, except that one needs now to introduce a replica matrix Q ra ,sb 
where the elements Q ra ,ra are fixed to the value qc for all r, a. In terms of Q and Ti the 
free-energy has the form fl32"|) : 



lim S.P.- 

n^O n 



f l V Q p 
4 ^ 



a,sb 9 



r,,,sbQra,sb + ~2^ 1 ^ 1 ^ 



(39) 



ra,sb ra,sb 

Supposing an ansatz such that Q\ a ,ib = Q2a,2b and Qi a ,2b = Q2a,ib, all to be taken as Parisi 
matrices, it is useful to define the submatrices Q a b = Qi a ,u and P a b = Qi a ,2b- 
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Before considering the explicit solution of fl39|), we just mention the difference in procedure 
to calculate the free energy /oT(<7cO- In this case, by the definition (p3|), the h fields are 
statistically independent for the two copies, and within each copy h a = h b if a and b belong 
to the same group of m replica, while h a and h b are independent in the other case. With 
the same meaning as before, we can then write: 



H, 



al) 



p/2 a=b 
p/2 g p_1 a,b in same group of m replica 
a,b in different groups 



(40) 



A 



ah 



Va,6. 



(41) 



keeping this difference in mind, the formal expression for the free energy f^fiqc) is still 
given by fl39|). 

We now want to find solutions for the matrix Q ra s b in both cases. 



A. Analysis of the saddle point equations 

We would like now to show the existence of solutions to the saddle point equations that 
have respectively V[f(qc) > for qc ^ q and V^ p (qc) > for qc ^ 0. 

The simplest Ansatz that we can try for the matrices Q a ^ and P a ^ is of 1RSB type. 

g = < x < mi 
q(x) = { (42) 

<?i mi < x < 1 

p(x)= <Po = 0<x< mi m 
Pi rrii < x < 1 

As we will see below, this is a valid Ansatz giving a positive potential for the case of V^^gc)- 
Unfortunately, we numerically found that the best variational solution of this form for 
fiTilc) gives negative values of the potential in some range of qc- A good positive solution 
is obtained at the level of 2RSB, which is the form we use, solving numerically the variational 
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problem with respect to the parameter mi, q±, q2,Pi,P2 of the functions: 



q(x) 



p(x) 



q = < x < m/2 
qi m/2 < x < mi 
q2 TXli < x < 1 

Po = < x < m/2 
Pi m/2 < x < mi 
p 2 mi < x < 1 



(44) 



(45) 



that parameterize respectively the matrices Q and P. With both parametrizations ([42"l , f43D 
and (|4|,|45]), the saddle points in eq. ([59|) turn out to be maxima with respect to all the 
variational parameters. The resulting expression for the potentials in terms of the variational 
parameters can be easily obtained substituting (f0],[I^) in fl39|); given its length we do not 
reproduce it here. 



B. Numerical evaluation of the variational potentials 

Before analyzing the potentials for generic values of t it is useful to discuss its behavior 
for t — 1 |L3], |l|, [L6|], where V^i and Vo,i coincide. 

At temperatures higher than the "dynamical transition temperature" |23| T^, a single, 
locally quadratic absolute minimum in q = is present, with V^f (0) = V^f (0) = 0. When 
the temperature is lowered below Td but is still higher than the thermodynamical transition 
temperature T s , besides the q = 0, V\ r f (0) = 0, a second local minimum, higher than the 
previous one appears for a value of q > 0. When T reaches T s and for values of T < T s , 
the replica symmetry is broken, the two minima are degenerate V^^O) = V^{q) = and 
the argument of the high overlap minimum q equals to the solution of the 1RSB variational 
problem at that temperature. 

The behavior of the t < 1 replica potentials is investigated solving numerically the vari- 
ational equations with respect to all the variational parameters. 

Let us first discuss the simple case of the high temperature replica symmetric phase for 
T > T s . In this case only the V^ p potential has to be considered and we show in fig. |l| that 
the required inequality V^ p (gc*) > for qc ^ 0, already satisfied for t — 1, remains valid 
for every t < 1. 
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Let us then analyze the RSB phase (T < T s ) when t < 1 and show that the minimum in 
zero (case 1) or in q (case 0) becomes respectively higher than the absolute one. 

We show in fig. [5] the results of the saddle point evaluation of the potential V^ p (gc) f° r 
different values of the parameter t at a fixed temperature T < T s , having verified that the 
behavior does not change for different values of temperature below the static transition. 

In fig. ^|we plot the potential V[l P {qc) for different values of the parameter t at the same 
temperature T < T s as in fig. |[ 

The two potentials V™ p (qc) and V^ ep {qc) are greater than zero for every value of qc 
different, respectively, from q and zero, and are quadratic near the (lower) minimum. 

We numerically find a quadratic dependence on 1 — t, for t ~ 1, of the value of the 
potential in the higher local minimum (fig. [| ^), implying that the constant 1/C(t) behaves 
as (1 — t) 2 for t close to one. 

V. LOWER BOUND FOR THE INTERPOLATING FREE ENERGY. 

To conclude the argument we now show that the constrained free-energies calculated in 
the previous section within the replica approach indeed supply lower bounds for the true 
ones. 

We just repeat here the argument given by Talagrand, specializing to our simple case. For 
each value of t we define an interpolating Hamiltonian H u that relates, varying the parameter 
u in the interval [0, 1], the two replica interpolating system with a suitable paramagnet. 

Once again, we only discuss the case of fij{qc), the same procedure being valid with 
minor modifications for fo t t{qc)- 

We have found in the previous section a good positive solution at the level of 2RSB; to 
deal with this case within an interpolating scheme we need three families of fields for each 
copy. One then defines: 

N 

H u [a, r] = VtdH p [a] + VtdH p [r] - y/1 - tj^ (( h i + h \i)°i + (hi + h%dn) 

N 

- y/(i- u)t (to!,* + 9l,i + glM + (g 2 i,i + gl,i + gidn) , (46) 
i=i 

where the fields h and hd are defined as in (||) and h\ and h\ chosen independently, while 
the variances and the covariances of the fields g can be written, in a notation consistent 
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with section ITVL as: 



P p-i 



P, 



l-q P 2 



1,2 



E 9 [gig 2 ] = \p\ 



E g [9l9l]=lf,- 1 E 9 [g\ g ^=\W^-f 2 



In a completely analogous way as in (||), we define the free energy: 

2 



fu 



(3Nm 



Ej\n 



Eg\,glh [E g i, g %Zc,u 



(47) 



(48) 



Z, 



C,u 



VaVTE h i :h 2 g i g 25(q(cT, r) - q c ) exp -pH u [a, 



where it is evident that f u =i,t = ft(lc)- On the other hand, one can check that f u =o,t 
reproduces the non trivial part of the replica expression of ft{qc) and one has 

ft ep (qc) = fu=o,t+ t ^ P 2 - {l - (1 - m x )ql - (mi - ^-)q p + q p c - (1 - m^ - (mi - 

(49) 

when the parameters that appear in the variances (f47|) maximize the expression ( j49|) for 
t = 1. 

We can thus repeat the argument as for the "standard" interpolating method: by the 
usual integrations by part one sees that the derivative can be written, for every t, as: 

djv_ = t/3(p 1) m ^jp _ ^ _ + gg, — (1 — mi )pg - (mi - y )K} + ^« , 

(50) 



where 1Z U is a remainder with the explicit expression: 



K u = y(l-mi)<2g((T,r) p -pg((T,r)(gf 1 + p p 2 1 



9i 



p— l p— l 



fx ) + (p-^ilZ + Pl- £ti-fi)) 2 u 



+ 



t/3 



2 -(mi - —) (2q(a, rf - pq(a, r)^ 1 + ff 1 ) + (p - l)(gf + pj))^ + ( ? (<r, r)"> J51) 

The inequality f u= o = flTiQc) < fx,t(Qc) — fu=x follows from the positivity of 1Z U , that is 
evident from the convexity of the function q p for even p, in the same way as in (0). In the 
Appendix we give the expressions for the averages that appear in the remainder. 



VI. CONCLUSION. 



The main point of this paper is to present the Talagrand proof of the Parisi anzatz in 
disordered mean field models in the physicist's language. With the aim of emphasizing 
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the conceptual aspects of the proof, we have specialized it to the case of the spherical p- 
spin model, where the simplicity of the Parisi solution allows a great reduction of technical 
difficulties. The main advantages of our approach is that we circumvent the most involved 
part of the demonstration by a direct solution of a variational problem, which in the case 
faced in this paper allows for an easy numerical solution. In doing that, we emphasize the 
connections of the Talagrand proof with the well known method of coupled replicas and the 



construction of glassy potential functions previously known in the physical literature |L2 



Though our paper does not really contain new mathematical results, we hope that it could 
be useful to spread the Talagrand's results in a larger community. 

One of the difficulties of the Talagrand paper is that the proof was provided in a general 
Parisi scheme with arbitrary steps of RSB. This is needed in the case of the SK model or in 
the low temperature phase of the Ising p-spin model. There, the derivation is complicated by 
the need to go to the limit of continuous replica symmetry breaking. In particular, Talagrand 
showed that, for every t < 1, it is possible to consider a number k = k(t) of steps of RSB 
such that an interpolating system with j > k steps is solved by the Parisi Ansatz with j 
steps. One then has to show the positivity of j different potential functions. An alternative 
is to work directly in the continuous RSB limit, and to look for a lower bound for a suitable 
family of potential functions depending on a continuous index. These potentials can be 
computed explicitly in spherical models with full RSB |Q, and we hope to report soon on 
that case. 

Thanks to the more physical approach followed in this revisitation of the Talagrand proof, 
it can also be easier to analyze more involved mean field systems, like those diluted, where a 
lower bound for the free energy was recently found || § by an extension of the interpolating 
procedure. 
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VII. APPENDIX. 



Generalizing the procedure used to define the averages fllTD , the first step to define (-) 0u , 
(■) lu and (-) 2u is to consider three four-replica Hamiltonian with the general form: 



N 



H = VutH p [a] + VutH p [r] - y/l + h \i)°i + + h d,i) T i) 

N 

- y/(l-u)t^2 ({gl ti + + g l d M + (<?i + <??2 + 4<) r 

i=l 

JV 

+ v^t^M + - (& + + & + ~ h dM] 

i=l 
N 

- - u)t {{g{,i + 92,i + 9dM + (ati + 9I2 + 9h)r-) , (52) 

8=1 

where the new introduced fields have the same distribution of the corresponding h and g 
fields: h have the same distribution of h, and g the same distribution of g. 
The three different Hamiltonians are specified by: 

1. H : the h and g fields are respectively independent of the fields h and g. 

2. Hi: g\ = g\, g\ = g\ and h = h, while the other h and g fields are independent of the 
corresponding h and g fields. 

3. H 2 : g\ = g{, g\ = gf, h = h, g\ = g\, ~g\ = g\\ the "annealed" fields h r id , h r id , g\ d and 
9 id ( r = I; 2) are always independent. 



We can now define, as done in (M), the Boltzmann averages of a given function of two spin 
configurations k(a, r), averaged over all the annealed fields h r id , h r id , g r id , g"!- d , r = 1,2: 

n(k(a, r)) 2 ,i )0 = " d " d " d _ r = 1, 2 , (53) 

^2,1,0 

/+OO 
^r, A or, r VaVTS(q(a, a') - q c )S{q{T, r') - g c ) exp (-/3H 2 ,i,o) • 
1 ,d i,d^i,d ,:j i,d 
■OO 
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The three averages that appear in the remainder are then defined by: 



E, 



(k(<T,T))„ = Ej- 



(E 9W Z?) 2mi " E gW2 Z 2 ^n(k(a, r)) 2 



E n 



E n i n 2Zo~ 



7,1 o 2 h I-Jn 1 n 2 -£'2 



(k(a,r)) hu = Ej- 



E 9\ 



/ \ 

I 771 -7 mi \ 4m l 



(54) 



(55) 



(^ 9 i, 92 2 ,^,9 2 2 ^o mi ) ^i, 9 2 i§ i i§ 2Z mi fi(A;(cT, r)) 



E 9i ,g\ ,h,9l,glh{ E g\ ,g\ ,~g\ ,g 2 2 Z o 



mi 



(56) 
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Figure captions 



Figure 1: The potential V^ p (qc) for different values of t at the temperature T s = 0.503 < 
T = 0.526 <T d = 0.544, for p = 4. From the bottom to the top the parameter t runs from 
1 to 0, decreasing by 0.1 at each plot. We have plotted only the results for q c > 0, being 
the potential symmetric for negative values of qc- 

Figure 2: The potential V[f(qc) at the inverse temperature /3 — 3 and p — 4. From the 
bottom to the top the parameter t runs from 1 to 0.9, decreasing by 0.01 at each plot. In 
the zoom is shown the potential near the higher minimum at the same inverse temperature 
and for the same values of the parameter t. Note that the potential is not symmetric with 
respect to positive and negative values of qc- 

Figure 3: The potential V^ p (?c) at the inverse temperature (3 = 3 and p = 4. From the 
bottom to the top the parameter t runs from 1 to 0, decreasing by 0.1 at each plot. In the 
zoom is shown the potential at the same inverse temperature near the higher minimum, for 
t varying from 1 (at the bottom) to 0.9 (at the top) decreasing by 0.01 at each plot. Note 
that the potential is now symmetric with respect to positive and negative values of qc, being 
the minimum at g i2 = 0. 

Figure 4: Value of the potential V[f in the higher local minimum, as a function of 1 — t, 
at the inverse temperature /3 — 3 and p — 4. The continuous line is a quadratic fit with the 
function f(x) = ax2, a = 2.26684. 

Figure 5: Value of the potential V^ p in the higher local minimum, as a function of 1 — t, 
at the inverse temperature (5 — 3 and p — 4. The continuous line is a quadratic fit with the 
function f(x) = ax 2 , a = 0.149782. 
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